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Abstract 

An approach is proposed to obtain some exact explicit solutions in terms of the 
Weierstrass' elliptic function p to a generalized Benjamin-Bona-Mahony (BBM) 
equation. Conditions for periodic and solitary wave like solutions can be expressed 
compactly in terms of the invariants of p. The approach unifies recently established 
ad-hoc methods to a certain extent. Evaluation of a balancing principle simplifies 
the application of this approach. 
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1 Introduction 

The Benjamin-Bona-Mahony equation has been investigated Kularized 
version of the Korteweg-de Vries equation for shallow water waves [1] . It incor- 
porates nonlinear dispersive and dissipative effects [2,3]. In certain theoretical 
investigations the equation is superior as a model for long waves, and the 
word "regularized" refers to the fact that, from the standpoint of existence 
and stability, the equation offers considerable technical advantages over the 
Korteweg-de Vries equation. In addition to shallow water waves, the equa- 
tion is applicable to the study of drift waves in plasma or the Rossby waves 
in rotating fluids. Under certain conditions, it also provides a model of one- 
dimensional transmitted waves. These find applications in semiconductor de- 
vices, optical devices, etc. [1]. - Apart from these applications solutions to the 
BBM equation are interesting in and of themselves. 

Recently, Wazwaz [4] has introduced five new ansaetze, sinh-cosh-ansatz I - 
III, a sine-cosine- ansatz and a tanh-sech-method to find exact solutions to a 
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generalized BBM equation 

ih + u x + a u n u x + u xxx = 0, n > 1, (1) 

with constant parameter a. The ansaetze are direct approaches and "ad hoc" 
as a method, but yield (according to [4]) new results. Therefore, it seems rea- 
sonable to introduce and apply a further direct method to the above equation. 
As an interesting "feature" this method has compact conditions for periodic 
and solitary solutions. The power of the method is its ease of use and the 
existence of the mentioned conditions. Furthermore, it will be shown that the 
sinh-cosh-ansatz I-II, the sine-cosine-ansatz, the tanh-sech-method and the 
sinh-cosh-ansatz III for n — 1 are subcases of the method presented here. 



The main features of this method are outlined in sec. 2 and applied to the BBM 
equation in brief. A balancing principle is evaluated in subsec. 2.1. Periodic and 
solitary solutions to the BBM are presented in subsec. 2.2. The relationship 
between this method and the ansaetze introduced in [4] are shown in sec. 3. 
A summary is given in sec. 4. 



2 Outline of the method [5,6] 

The starting point is a nonlinear wave and evolution equation (NLWEE) 

rKx) )f i'(x) ) «"(x),... 1 /)(x) r .. lI1 W( x )| = { W) (2) 



where r is a function of u(x) and its partial derivatives, the independent vari- 
able x has n components Xj, u^(x) denotes the collection of mixed derivative 
terms of order k and t(u) is a trigonometric function in or t(u) = 0. For no- 
tational simplicity the independent variables x, t will be used in the following. 
Equation (2) describes a certain dynamical system by means of a (wave) func- 
tion m(x). An travelling wave ansatz 

M 

u(x, t) - g(f(z)) = £ ai f(zY, MeN, (3) 

i=0 



where / is supposed to obey the nonlinear differential equation 

f^r) 2 = af + 4/3/3 + 67/2 + 46 f + e = R{f) ' (4) 



2 



(with real a, (3, 7, 5, e, z — fi{x — ct), f(z)) transforms the NLWEE into 
an polynomial equation P(f) = 0. Vanishing coefficients in the polynomial 
equation P(f) = lead to equations which partly determine the coefficients 
a, (3, 7, 5, e in eq. (4). In general, the coefficients depend on the structure and 
parameters of the NLWEE and, finally, on the parameters of the transforma- 
tion ip — > g. The parameter M in eq. (3) can be derived by balancing ( cf. 
subsec. 2.1). 

Thus, the problem of finding a solution to the NLWEE is reduced to finding 
an appropriate transformation that leads to eq. (4), which, in this sense, is the 
basis of the following analysis (for reference purposes it is called the "basic 
equation" of the associated NLWEE). 



As is well known [7,8] the solutions f(z) of (4) can be expressed in terms of 
Weierstrass' elliptic function p. It reads 

/( " ) = /0 + 4b(,^J 3 )- ) ^U)]' (5) 



where the primes denote differentiation with respect to / and /o is a simple 
root of R(f). 

The invariants g 2 , #3 of Weierstrass' elliptic function p{z\ g 2 , #3) are related to 
the coefficients of R(f) by [9] 

92 = at - A(35 + 3 7 2 , (6) 

g 3 = a 7 e + 2f3~f5 - a5 2 - 7 3 - e(3 2 . (7) 

The invariants and the discriminant (of p and R [9]) 

A = g\ - 27gl (8) 

are suitable to classify the behaviour of f(z). The conditions [11] 

A^O or A = , 92 >0, 53 >0. (9) 

lead to periodic solutions, whereas the conditions 

A = 0, g 2 >0 , <0 (10) 

are associated with solitary solutions. 
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If A = 0, the solution (5) can be simplified because <p(z;g2 1 g 3 ) degenerates 
into trigonometric or hyperbolic functions [13, pp. 651-652]. 

Thus, in this case, periodic solutions according to eq. (5) are determined by 

R'(fo) 



M = f + 



-f -^ + |e lC sc^,) 



, A = 0, g 3 > 0, (11) 



and solitary wave like solutions by 



f(z) = f + 



f(z) = f + 



R'(fo) 



24 - R"(f )zi 



, A = 0, g 3 < 0, (12) 
, A = 0,^3 = 0, R"(fo)<0, (13) 



d - ™^ + 3eicsch 2 (v^^) 
6R'(fo)z 2 



where e± = \ \g 3 \ in eq. (11) and e\ 



1 3 



g 3 \ in eq. (12). 



In general, f(z) (according to eq. (5)) is neither real nor bounded. Conditions 
for real and bounded solutions f(z) can be obtained by considering the phase 
diagram of R(f) [12, p. 15]. They are denoted as "phase diagram conditions" 
(PDC) in the following. Examples of a phase diagram analysis is given in refs. 
[10,11,12]. 



2. 1 Balancing principle 



As mentioned above, the exponent M in eq. (3) is determined by use of a 
balancing principle. In this subsection a balancing principle for the approach 
outlined above is evaluated. 

As a consequence of the transformation u — > g(f(z)) (cf. eq. (3)), the NLWEE 
can be rewritten in terms of g(f) and its derivatives. To determine the leading 
M exponent 



g(f) = «o + a M f 



M 



(14) 



has to be substituted into the NLWEE in question that leads to an expression 
of the form 



P(f)=0, 



(15) 
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VW 4 + W 3 + 6 7 / 2 + W + e g(/) = o, 
+ vV 4 + W 3 + 6 7 / 2 + 45/ + e Q(f) = 0, 



(16) 
(17) 



where P and Q are polynomials in /. To determine M in ansatz (3) - according 
to type (15), (16), (17) - the polynomial equations P(f) = 0, Q(f) = or the 
equation 



P(f) 2 - R(f) Q(f) 2 = o 



(18) 



is considered and every two possible highest exponents are equated. To avoid 
applying this procedure for every NLWEE anew, the highest possible expo- 
nents of the polynomials that occur for special derivatives in the nonlinear 
wave and evolution equation in question have been evaluated. One obtains 
from equation (4) 



d 



d 



_ _> j a p + A(3 p + 67/ 2 + 4Sf + 
dz 3 



d 



(19) 
(20) 



(6a/ 2 + 12(3 f + 6- )\ n/' • 1 //"• • Cr P ■ U)f ■ <^ (21 ) 



+ (2a/ 3 + 6(3 f + 67/ + 45) yW 4 + 4/3/ 3 + 6 7 / 2 + 45/ + e^. 



Consideration of these relations and eq. (14) leads to the following list of the 
highest possible exponents (HPEs) (cf. tab. 1). 



Function 


9 


9z 


3^ 


9zzz 


3"3 2 




HPE of P 


M 





M + 2 










HPE of Q 





M - 1 





M + l 


M (n + 1) - 1 





Table 1 

Highest possible exponent (HPE) of function g and its derivatives according to (14) 
and equations (19) - (21). 



2.2 Solutions to the BBM equation 



Using a transformation 1 

1 Balancing g zzz and g n g z yields M = | according to tab. 1. To simplify the 
comparison with the ansaetze in [4] M = ^ is chosen. It is obvious, that solutions 
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u(x, t) = g(z), z = fj,(x — ci), 
g(z)=h(z)$, (22) 

2 
3=0 

where / is supposed to obey eq. (4), leads to a polynomial equation P(f) = 0. 
Setting the coefficients equal to zero yields 



ao2n 2 aain 2 



7 : 



4(n + 4)(n+ l)^ 2 ' ^ 8(n + 4)(n + l)^ 2 ' 
n 2 (-3a(a 2 + 4a a 2 ) + a 2 (c - l)(n + 4)(n + 1)) 
96a 2 (n + 4)(n+ 1) ; 



aX^a 2 - 12a a 2 ) + a 2 (c - l)(n + 4)(w + 1)) 

64a 2 (n + l)(n + 4)^ 2 ' 1 j 

n 2 (a(a 4 - 12a a 2 a 2 + 48aga 2 ) + a 2 (a 2 - 8a a 2 )(c - \){n + 4)(n + 1)) 



e — 



a = — V a = 

4a 2 4aa 2 



2 



64a|(n + l)(n + 4)^ : 
aa 2 + a 2 (c-l)(n + 4)(n+l) 



According to eqs. (6) - (8) the invariants of the Weierstrass' elliptic function 
and the discriminant read 



C-l V C-l 3 n 6 . r\ • r <A 

92 = - ^-,03 = -- —z; A = Oif a = — ; (24) 

y2 3072/x 4 ' y3 884736^' 4a 2 ' V ; 



_ (c-l) 2 n 4 (c-l) 3 n 6 

92 ~ 192^ ,93 ~ 13824^ 6 ' 

aa 2 + a 2 (c-l)(n + 4)(n+l) 

ito = . (25) 

4aao 



As can be seen from eqs. (24) and (25), periodic solutions are given if c < 1, 
whereas c > 1 leads to solitary solutions, according to conditions (9), (10). 

2 

If a = periodic and solitary solutions can be evaluated according to eqs. 



(26), (27) and (28) can easily be rewritten as h(z)n as expected from the balancing 
result. 
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(22), (11) and (12) 2 , respectively. The periodic solution reads 



( C -l)(n + 4)(n + l)8ec 2 [g(x-rf)^/(l- C )] 
u(x, t) — <i ^ f . (2b) 



Solitary solutions are given by 



(c - l)(n + 4)(n + l)sech 2 [f (x - ct)y/(c - 1)] \ " 
u(x, *) = <( — 1 } , (27) 



'{aa? + 2a 2 ^2aa 2 (c-l)(n + 2)(n+l)sech[f(x - ct)y/{c - 1)]} 2 \ " 
u( *'' )= l 16^ - (28) 



If a = aa i+ a2 ( c 4 ^" +4 ' )( ' n+1 ' ) the basic equation has two double roots. Therefore, 
the general solution to eq. (4) [7,8] 

f(\ = f, \f^M M2 § Ml + ^U)b(^2,<? 3 ) - ±R!'(fo)] + ±RUo)R"\h) 

where Jo is an y constant, not necessarily a simple root of R(f), has to be 
considered. Evaluation yields exactly solutions (26) and (27). 



Equations (26), (27) and (28) present periodic and solitary solutions to the 
BBM equation. Because of the generalized ansatz there are no restrictions for 
a\ and a-i- This may be advantageous to adapt the theoretical results to a 
physical problem that is modelled by the BBM equation. 



3 Relationship between the approach above to some other ad-hoc 
methods 



Recently, Wazwaz has introduced five ansaetze to solve the BBM equation [4]. 
Comparing these ansaetze with the approach outlined above, shows, that four 
of the ansaetze occur as special cases of the approach here. The sinh — cosh- 
ansatz III can only be classified as a special case if n — 1. It is remarkable, 
that this classification as "special cases" does not depend on the nonlinear 

2 The case #2 = 93 = is neglected here, because according to the PDC there exist 
no bounded solutions. 
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equation in question, but is a characteristic of the approach outlined above. 
Therefore, features like the balancing principle, conditions for periodic and 
solitary wave solutions and the PDC can also be applied to the ansaetze given 
in [4]. This may be of interest because these ansaetze have also been applied 
succesfully to the Boussinesq equation [14]. 

As mentioned above the exponent M in eq. (3) can always be obtained by 
balancing 3 . Therefore, the base of the ansaetze in [4] has to be compared 
with eq. (4). The comparison is as follows: An ansatz f(z) according to [4] 
with z = x — ct is considered. The function / is differentiated with respect 
to z and (J^r) * s rewritten in terms of /. Comparing this expression with 
eq. (4) leads to the relationship between the parameters in the ansatz and the 
parameters in the basic equation. The results are presented in tab. 2. 

It should be pointed out that the following classification of the ansaetze as 
subcases of a generalized ansatz (3), (4) has several advantages: Firstly, it 
is shown that ad-hoc methods can be unified to a certain extent. Secondly, 
periodic and solitary solutions can be deduced systematically according to 
conditions (9) and (10). Thirdly, certain procedures, e. g. balancing, has only 
to be done once and not for every ansatz anew. Finally, some of the remarks 
concerning solutions evaluated in ref. [4] (confered to as "W") can be approved 
from another point of view: The remark, that solutions given by eqs. (W20), 
(W21) are consistent with solutions (W18), (W19) can be approved by consid- 
ering tab. 2. For A = ±1 in the cosh I -ansatz and A = ±i in the sinh I-ansatz 
the parameters of the basic equation are equal (a = in both cases), so that 
the same solutions have to be expected. 



4 Summary 

A method is proposed to obtain exact elliptic solutions to NLWEEs. The 
method is a generalization of some recently established ad-hoc methods [4] (cf. 
tab. 2) and, furthermore, includes conditions for periodic and solitary solu- 
tions. These are given compactly in terms of the invariants of the Weierstrass' 
elliptic function p. Periodic and solitary solutions can be deduced system- 
atically from the general solution (5) and (29), respectively. Thus, a further 
promising ad-hoc method has been shown and, exemplary, solutions to the 
BBM equation have been evaluated. A balancing principle has been evaluated 
that simplifies the application of this method (cf. tab. 1). It seems that this 

3 In ref. [4] the exponent in the sinh — cosh-ansaetze I-III is directly given as ^. 
This result is obtained by balancing here. 
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Name 


Ansatz f(z) 


Parameters of the basic equation 


Comment 


cosh I 


b 

l+A cosh(/i z) 


a=^(l-A 2 ),/3 = -g, 7 = ^ = e = 


special case of eq. (4) 


sinh I 


b 

l+A sinh(/i. z) 


a=^(l + A 2 ),/? = -g, 7 = ^,5 = e = 


special case of eq. (4) 


cosh II 


b 

l+A cosh 2 (fi z) 


a = + A), = -4(2 + A), 7 = ¥> 5 = e = 


special case of eq. (4) 


sinh II 


b 

l+A sinh 2 (p z) 


a = ^(1 - A),/? = -4(2 - A), 7 = ¥^ = e = 


special case of eq. (4) 


cosh III 
n = 1 


b cosh 2 (/j, z) 
l+A cosh 2 (p z) 


4AV(1+A) fl M 2 (36A 2 +26A) 
7 _ 2 ,WA+ b 2 ) )5 _ b ^ e _ Q 


special case of eq. (4) 


sinh III 
n — 1 


b sinh 2 (p z) 
l+A sinh 2 (p z) 


4AV(1-A) fl ^ 2 (-36A 2 +26A) 
« - b 2 ,fJ— b 2 

7" V( 1? +62) ,5-V, e -0 


special case of eq. (4) 


cosine 


A cos(/iz) 


a = /? = < J = 0,7=-^,e = AV 


special case of eq. (4) 


sine 


A sin(yuz) 


a = /? = < J = 0,7=-^,e = AV 


special case of eq. (4) 


tanh-sech 


tanh(/iz) 


a = ^ 2 , /? = 5 = 0, 7 = — 4) e = /i 2 


special case of eq. (4) 



Table 2: Relationship between the approach outlined here 
and some ad- hoc methods given in [4]. 



method is a useful and easily manageable tool with interesting features (bal- 
ancing principle, conditions for periodic and solitary wave solutions, PDC) to 
find out exact solutions to nonlinear differential equations, especially solitary 
ones. It may be advantageous that this quite general method can lead to free 
parameters in the solution, as shown for the BBM equation. 
It has recently been shown [5] for the Kadomtsev-Petviashvili equation and 
the nonlinear Schrodinger equation that if A 7^ and special conditions hold, 
eq. (5) can be rewritten in terms of Jacobian elliptic functions and used as 
a start solution for a linear superposition principle that enlarges the solution 
set of periodic solutions [15]. For the Novikov-Veselov equation it has been 
shown previously [16], that a 2-solitary wave solution can be evaluated by 
linear superposition of two 1-solitary wave solutions to the Korteweg-de Vries 
equation. These 1-solitary wave solutions have been evaluated by the approach 
outlined above. Thus, this approach can also be used as a basis to construct 
periodic superposition solutions and multi-solitary wave solutions. 
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